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O I Abstract 

Within the Hamiltonian approach in Coulomb gauge the ghost and gluon propa- 
gators are determined from a variational solution of the Yang-Mills Schrodinger equa- 
tion showing both gluon and heavy quark confinement. The continuum results are in 
^^ I good agreement with lattice data. The ghost form factor is identified as the dielectric 

function of the Yang-Mills vacuum and a connection between the Gribov-Zwanziger 
scenario and the dual Meissner effect is established. The topological susceptibility is 
calculated. 

Oh: 
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1 Introduction 
> 

Cf^ , The aim of the talk is the microscopic description of infrared properties of QCD like confine- 
2P I ment. We would like to see, for example, the emergence of the colour flux string between 
static colour charges. For this purpose, I will use the Hamiltonian approach to Yang- 
Q ■ Mills theory in Coulomb gauge. The organisation of my talk is as follows: In section 2, I 
00 , will briefly summarise the essential ingredients of the Hamiltonian approach to Yang-Mills 



O 
> 



theory in Coulomb gauge. Then I will present a variational solution of the Yang-Mills 
Schrodinger equation in section 3, which will result in a set of coupled Dyson-Schwinger 



/\ I equations. I will present the analytic solutions to these equations in the infrared and ultravi- 
c^ ■ olet and the numerical solution for the full momentum range. The resulting propagators will 
then be compared to the available lattice data. Then I will focus on two non-perturbative 
properties of the Yang-Mills vacuum: the dielectric constant in section 4 and the topological 
susceptibility in section 5. Finally, a summary is provided. 

2 Canonical quantisation of Yang- Mills theory 

In the canonical quantisation approach the gauge fields v4'^(x) are considered as the (carte- 
sian) coordinates and the corresponding conjugate momenta are defined by 
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where S is the action of the Yang-Mills field. The explicit calculation yields 

n^(x) = E^ix) , ns(x) = , (2) 

where E^{x) is the colour electric field. To avoid the problems arising from the vanishing 
temporal component of the canonical momentum, one imposes Weyl gauge ^o(^) ~ 0- '^^^ 
Yang-Mills Hamiltonian is then given by 

H = ]^j d?x{n\x) + B\x)) . (3) 

The canonical quantisation is carried out in the standard fashion by imposing the canonical 
commutation relation [A?(a;), n^(|/)] = 5ij5°'^5^{x — y), which promotes the canonical mo- 
mentum to the operator n"(x) = 6/i6A°;{x). By imposing Weyl gauge one loses Gauss' law 
from the Heisenberg equation of motion and Gauss' law has to be imposed as a constraint 
on the wave functional 

D ■ U{x)i!{x) = -gpm{x)i'{x) , (4) 

where pm{x) is the colour charge density of the matter fields and D^^ = 6"'^di + gA'}'^ 
^j^ab _ j-acbj^c^ jg ^^^^ covariaut derivative in the adjoint representation of the gauge field 
with f"-'''^ being the structure constant of the gauge group. The operator on the left hand side 
of Gauss' law is nothing but the generator of time-independent gauge transformations and 
in the absence of external colour charges, pm{x) = 0, Gauss' law expresses the invariance of 
the wave functional under space-dependent but time-independent gauge transformations. 

Instead of working with explicit gauge invariant wave functionals it is more convenient to 
explicitly resolve Gauss' law by fixing the gauge. Coulomb gauge is a particular convenient 
choice for this purpose. We implement the Coulomb gauge, d ■ A = 0, in the standard 
fashion into the scalar product of the wave functionals by means of the Faddeev-Popov 
method 

{^l^\0\cj>) = [ DA^JiA^)riA^)0[A^]<l>{A^) , (5) 



where 

J = Det(-£> ■ d) (6) 

is the Faddeev-Popov determinant. While in Coulomb gauge the gauge field is transversal 
the momentum operator 11 = II" + 11 contains both longitudinal II" and transversal 11 
parts. Resolving Gauss' law for the longitudinal part of the momentum operator yields 

n'l V = gd{-D ■dr'ptlj,p = pg + pm, (7) 

where 

p, = A^ ■ n^ (8) 

is the colour charge density of the gauge field. With this result the Hamiltonian in Coulomb 
gauge is found to be 

H = l [ (j-'n^jn^ + B^) + He , (9) 



where 

1 f .-1^,1 .^„ 9' 



He = -Jj-'n^^jn\\ = ^-Jj-'pi-D-d)-\-d'){-D-d)-'jp (10) 

is the so-called Coulomb Hamiltonian, which arises from the longitudinal part of the kinetic 
energy after resolving Gauss' law. The Hamiltonian (Q was first derived in Ref. [1]. 
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Figure 1: Ghost form factor d{k) (left) and gluon energy uj{k) from the variational solutions 
presented in [2]. 



3 Variational solution 

We wish to solve the Schrodinger equation 

H'iIj = EiIj 
by the variational principle 



\H\ 



rain 



:ii) 



(12) 



with suitable ansatze for the wave functional ip{A-^). This approach has been studied 
resently in Refs. [Sj H]. Inspired by the wave functional of a massless particle moving 
in a spherically symmetric potential in an s-state i}) = (p{r)/r, where r = J^''^(r) is the 
Jacobian of the transformation from the cartesian to the spherical coordinates for zero 
angular momentum we choose the following ansatz [3] 



^{A' 






(13) 



where the kernel uj{x,y) is determined from the variational principle flT2l) . In practice the 
so resulting equation for Lj{x,y) is converted into a set of Dyson- Schwinger equations for 
the gluon propagator 
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{At%x)Af\y)) = 6^%{x)-uj-\x,y) 



(14) 



did 



with tij{x) = 6ij — ^ being the transverse projector, and the ghost propagator 



G{x,y) 



D d 



-1 



{x\d{-A){-A)-'\ 



(15) 



Here we have introduced the ghost form factor d{—A), which describes the deviation of 
the QCD ghost propagator from the QED case, where d{—A) = 1. The resulting Dyson- 
Schwinger equations need renormalisation, which is well under control. Fig. [1] shows the 
solution of the Dyson- Schwinger equation for the gluon energy uj{k) and the ghost form 
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Figure 2: Left: Heavy quark potential given by eq. (TT8|) . Right: Running coupling constant. 



factor d{k), as shown in Ref. [2]. An analytic infrared and ultraviolet analysis of the 
Dyson-Schwinger equation shows the following asymptotic behaviour [H |5] 
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d{k) ~ k^ . 



(16) 



At large momenta the gluon behaves like a photon, which is in agreement with asymptotic 
freedom, while at small momenta the gluon energy diverges, which implies the absence 
of gluon states in the physical spectrum. This is nothing but a manifestation of gluon 
confinement. The infrared divergence of the ghost form factor is a consequence of the 
horizon condition 

d-\k = 0)=0, (17) 

which has been used as input in the renormalisation of the ghost Dyson-Schwinger equation. 
This is a necessary condition for the Gribov-Zwanziger confinement scenario. In fact, one 
can show there is a sum rule relating the infrared exponents of the ghost and the gluon 
propagator and an infrared divergent gluon energy requires also an infrared divergent ghost 
form factor, i.e. the horizon condition flTTl) . see Ref. 0. Fig. [2] shows the non-Abelian 
Coulomb potential 



V{\x-y\) = 9'({x\{-D . d)-\-d'){-D ■ d) 



<yc\x-y\ 



which for large distance indeed increases linearly [2] as the infrared analysis reveals. The 
Coulomb string tension ac sets the scale of our approach. Also shown in Fig. [2] is the 
running coupling constant which is infrared finite, for details see Ref. [5]. Fig. [3] shows the 
continuum results for the gluon energy and the ghost form factor in D = 2 + 1 dimensions 
[6] together with the corresponding lattice results, Ref. [7]. The agreement is not perfect 
but, given the approximation involved, quite satisfactory. 

In I? = 3 + 1 dimensions, previous lattice calculations performed in Coulomb gauge in 
Ref. [Bl E] showed an anomalous UV behaviour of the gluon propagator — YR : uj{k) r^ 
k^ , UV : uj{k) ~ fc^/^ — which is in strong conflict with the continuum result. However, 
one should mention that these lattice calculations assumed multiplicative renormalisability 
of the 4-dimensional gluon propagator, which give rise to scaling violations in the static 
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Figure 3: Results for the ghost form factor d{k) (left) and the gluon energy uj{k) in 2 + 1 
dimensions, as shown in [6]. 

propagator. Furthermore, these calculations did not fix the gauge completely, i.e. the 
residual time-dependent gauge invariance left after Coulomb gauge fixing was left unfixed. 
Furthermore, the Coulomb gauge fixing was done on a single time-slice, which is sufficient for 
the calculation of static (time- independent) propagators. However, one should keep in mind 
that in Coulomb gauge topologically non-trivial gauge configurations, which presumably are 
responsible for confinement, are discontinuous in time [lOj and as a consequence on a small 
lattice different results are obtained from different time slices. 

Recently, we have done improved lattice calculations with a complete gauge fixing [11]. 
In these studies, the energy dependence of the 4-dimensional gluon propagator could be 
explicitly extracted and it was found that the static gluon propagator is multiplicatively 
renormalisable and shows a perfect scaling. Fig. S] (left panel) shows the results for the 
gluon propagator of these calculations together with the continuum results. It is assumed 
here that the Coulomb string tension ac is identical to the string tension a from the Wilson 
loop. There is a very good agreement, in particular the ultraviolet and infrared behaviour 
matches perfect for lattice and continuum. What is also remarkable that the lattice result 
can be very well fitted by Gribov's original formula for the gluon energy 



uj{k) 




(19) 



with M = 0.88(1) GeV. 



4 The colour dielectric constant 



Consider the electric field generated by a charge density p in electrodynamics 

E = -a0,0=(-A)-V. 



(20) 



The longitudinal electric field resulting from the resolution of Gauss' law in the Yang-Mills 
case is given by a similar expression 
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Figure 4: Left: Lattice data for uj{k), compared to the solution of the Dyson-Schwinger 
equations. Right: Dielectric function e{k). 
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p = rf(-A)(-A)-V 



(21) 



except that the Green's function of the Laplacian is replaced by the ghost propagator ( jTSl) . 
The last expression has the form of the scalar potential in the presence of a dia-electric 
medium 



-A)-V 



(22) 



and the inverse of the ghost form factor d{k) can thus be identified as the dielectric function 
of the Yang-Mills vacuum 

e{k) = d-\k) . (23) 

Fig. m (right panel) shows the so defined dielectric function. It satisfies < e{k) < 1, which 
is a manifestation of anti-screening while in QED we have e{k) > 1, which corresponds to 
ordinary Debye screening. Furthermore, at zero momentum the dielectric function vanishes, 
showing that in the infrared the Yang- Mills vacuum behaves like a perfect colour dia-electric 
medium. The vanishing of the dielectric function in the infrared is not an artifact of our 
solutions of the Dyson-Schwinger equations but is guaranteed by the horizon condition, 
which is a necessary condition for the Gribov-Zwanziger confinement scenario. A perfect 
colour dia-electric medium e = is nothing but a dual superconductor. (Here, "dual" 
refers to an interchange of electric and magnetic fields and charges.) Recall in an ordinary 
superconductor the magnetic permeability vanishes /i = 0. This shows that the Gribov- 
Zwanziger confinement scenario implies the dual Meissner effect [12j . 



5 Topological susceptibility 



As first shown by Adler [13] and Bell and Jackiw [H], the Ua{^) symmetry is anomalously 
broken which gives rise to an extra mass term to the rj', which by the Witten-Veneziano 
formula 
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Figure 5: Topological susceptibility x as a function of the ratio 0^0. 



is expressed by the topological susceptibility 

which is the correlation function of the topological charge density 

^2 



q{x) 
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(25) 



(26) 



Furthermore, in eq. (!24l) Nf denotes the number of flavours and F^ ~ 93 MeV is the 
pion decay constant, x vanishes in all orders of perturbation theory and is thus an ideal 
observable to test the non-perturbative content of our vacuum wave functional. In the 
Hamiltonian approach one finds the following expression for the topological susceptibility 

m 



Vx 




fB\x)\0)-2Y: 



|(n|/B-n|0)p 



Er, 



(27) 



Here In) denotes the exact excited states of the Yang- Mills Hamiltonian with energies En- 
These eigenstates are of course not known. We work out the matrix elements in eq. fl27|) to 
two-loop order. In this order only two and three quasi gluon states 



af(x)af(y)|0),af(x)af(y)af(^)|0) 



:28) 



contribute where our vacuum state is annihilated by the operators a°(x), i.e. a°(x)|0) = 0. 
The resulting expression for the topological susceptibility is ultraviolet divergent and needs 
renormalisation. For this aim we exploit the fact that x vanishes to all order perturbation 
theory and renormalise the expression (!27|) for x by subtracting each propagator by its 
perturbative expression. This renders x (1271) finite. Furthermore, since the momentum 
integrals in this expression are dominated by the infrared part we replace the coupling 
constant, which, in principle, should be the running one, by its infrared value. The results 
obtained in this way for the topological susceptibility are shown in Fig. [5] (right panel) as a 
function of the ratio (Tc/ct- Choosing ac = l.Scr which is the value favoured by the lattice 
calculation [S] we find with y/a = 440 MeV 

X = (240 MeV)^ . (29) 

This value is somewhat larger than the lattice prediction x = (200 — 230 MeV)'^. 



6 Summary and Conclusions 

I have presented a variational solution of the Yang-Mills Schrodinger equation in Coulomb 
gauge using a Gaussian type of ansatz for the vacuum wave functional. We find a gluon 
energy which is infrared divergent, which is a manifestation of gluon confinement. Further- 
more, we have found a static colour charge potential which at large distances rises linearly, 
as one expects for a confining theory. The propagators calculated within this approach are 
all in satisfactory agreement with the lattice data. I have then shown that the inverse of the 
ghost form factor can be interpreted as the colour dielectric function of the QCD vacuum. 
The horizon condition, a necessary condition for the Gribov-Zwanziger confinement sce- 
nario to work, implies that in the infrared the QCD vacuum is a perfect colour dia-electric 
medium, which is nothing but a dual superconductor. In this way the Gribov-Zwanziger 
confinement scenario implies the dual Meissner effect. Finally I have presented results for 
the topological susceptibility calculated in the Hamiltonian approach with our vacuum wave 
functional. For reasonable values of the Coulomb string tension we find results close to but 
somewhat larger than the lattice data. The results obtained so far in this approach are 
quite encouraging for further investigations. A natural next step would be the inclusion of 
dynamical quarks. 
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